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Abstract 



The basic functional equations for connected and 2-connnected graphs can be 
traced back to the statistical physicists Mayer and Husimi. They play an essential 
role in establishing rigorously the virial expansion for imperfect gases. We survey this 
approach and inspired by these equations, we investigate the problem of enumerating 
some classes of connected graphs all of whose blocks are contained in a given class 
B. Included are the species of Husimi graphs (B = "complete graphs"), cacti (B = 
"unoriented cycles"), and oriented cacti (B = "oriented cycles"). For each of these, we 
consider the question of their labelled or unlabellcd enumeration and of their molecular 
expansion, according (or not) to their block-size distributions. 



Les equations fonctionnclles de base pour les graphes connexes ct 2-connexes 
rcmontcnt aux physiciens Mayer et Husimi. Ces relations sont importantes pour 
ctablir rigoureusement le developpemcnt du viriel pour les gaz imparfaits. Nous 
presentons cette demarche et inspires par ces relations, nous examinons lc probleme du 
denombrement dc quclqucs classes de graphes connexes dont les blocs sont pris dans 
une famille B donnee. Cela inclut les graphes de Husimi (B = "graphes complets"), 
les cactus (B = " polygones" ) , et les cactus orientes (B = "cycles orientes"). Pour 
chacune de ces especes, on s'interesse au denombrement etiquete, ordinaire ou selon la 
distribution des tailles des blocs, au denombrement non-etiquete, et au developpement 
moleculaire. 

1 Introduction 

1.1 Combinatorial species and functional equations for connected graphs 
and blocks 

Informally, a combinatorial species is a class of labelled discrete structures which is closed 
under isomorphisms induced by relabelling along bijections. See Joyal ^31 and Bergeron, 
Labelle, Leroux [2j for an introduction to the theory of species. To each species F are 
associated a number of series expansions among which are the (exponential) generating 
function, F{x), for labelled enumeration, defined by 
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n>0 



(1) 



where \F[n] | denotes the number of F-structures on the set [n] = {1, 2, . . . , n}, the (ordi- 
nary) tilde generating function F(x), for unlabelled enumeration, defined by 

F(x) = £ F n x n , (2) 

n>0 

where F n denotes the number of isomorphism classes of F-structures of size n, the cycle 
index series Zf{x\,X2,x%, • • •), defined by 

Z F (x 1 ,x 2 ,x 3 ,---) = Y / -A E fixi? H x l x 2X?'-\ , (3) 

n>0 n ' \<reS„ J 

where S n denotes the group of permutations of [n], fixF[u] is the number of F-structures 
on [n] left fixed by a, and aj is the number of cycles of length j in a, and the molecular 
expansion of F, which is a description of the F-structures and a classification according 
to their stabilizers and will be discussed later. 

Combinatorial operations are defined on species: sum, product, (partitional) compo- 
sition, derivation, which correspond to the usual operations on the exponential generat- 
ing functions. And there are rules for computing the other associated series, involving 
plethysm. See [2] for more details. An equality F = G between species is a family of 
bijections between structures, 

au ■ F[U] -> G[U], 

where U ranges over all underlying (labelling) sets, which commutes under relabellings. 
Such an identity gives rise to an equality between all the series expansion associated to 
species. 




Figure 1: A simple graph g and its connected components 

For example, the fact that any simple graph on a set (of vertices) U is the disjoint 
union of connected simple graphs (see Figure ^) is expressed by the equation 

9 = E(C), (4) 

where G denotes the species of (simple) graphs, C, that of connected graphs, and E, the 
species of Sets (in French: Ensembles). There correspond the well-known relations for 
their exponential generating functions, 

Q{x) = exp(C(x)) (5) 



2 



and for their tilde generating functions, 



g{x) = z E (c(x),c(x 2 ),...) 

= ex pfg^(* fe )j. ( 6 ) 

Definitions. A outpoint (or articulation point) of a connected graph g is a vertex of g 
whose removal yields a disconnected graph. A connected graph is called 2-connected if 
it has no cutpoint. A block in a simple graph is a maximal 2-connected subgraph. The 
block- graph of a graph g is a new graph whose vertices are the blocks of g and whose edges 
correspond to blocks having a common cutpoint. The block- cutpoint tree of a connected 
graph g is a graph whose vertices are the blocks and the cutpoints of g and whose edges 
correspond to incidence relations in g. See Figure El 




a) b) c) 



Figure 2: a) A connected graph g, b) the block-graph of g, c) the block-cutpoint tree of g 

Now let B be a given species of 2-connected graphs. We denote by Cg the species of 
connected graphs all of whose blocks are in B, called Cg-graphs. 

Examples 1.1. Here are some examples for various choices of B: 

1. If B = B a , the class of all 2-connected graphs, then Cb = C, the species of (all) 
connected graphs. 

2. If B = K2, the class of "edges", then Cg = (2, the species of (unrooted, free) trees 
{& for French arbres). 
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3. If B = {P m ,m > 2}, where P m denotes the class of size-m polygons (by convention, 
P2 = K2), then = Ca, the species of cacti. A cactus can also be defined as a 
connected graph in which no edge lies in more than one cycle. FigureOl a), represents 
a typical cactus. 




a) b) 



Figure 3: a) a typical cactus, b) a typical oriented cactus 

4. If B = K3 = P3, the class of "triangles", then Cb = 5, the class of triangular cacti. 

5. If B = {K n ,n > 2}, the family of complete graphs, then Cb = Hu, the species of 
Husimi graphs, that is of connected graphs whose blocks are complete graphs. They 
were first (informally) introduced by Husimi in |12| . A Husimi graph is shown in 
Figure El b). See also Figure It can be easily shown that any Husimi graph is the 
block-graph of some connected graph. 

6. If B = {C n ,n > 2}, the family of oriented cycles, then Cb = Oc, the species of 
oriented cacti. Figure |21 b) shows a typical oriented cactus. These structures were 
introduced by C. Springer |2H] in 1996. Although directed graphs are involved here, 
the functional equations (J7J) and (|12[) given below are still valid. 

Remark. Cacti where first called Husimi trees. See for example [201 an d |29j . 

However this term received much criticism since they are not necessarily trees. Also, a 
careful reading of Husimi's article |12j shows that the graphs he has in mind and that he 
enumerates (see formula (|42[) below) are the Husimi graphs defined in item 5 above. The 
term cactus is now widely used. See Harary and Palmer ^Ql- Cacti appear regularly in 
the mathematical litterature, for example in the classification of base matroids [20], and 
in combinatorial optimization 

The following functional equation (see (JJJ)) is fairly well known. It can be found in 
various forms and with varying degrees of generality in [2] , JO] , ECU , UH , CH| , HH , HU , 
[2*Sj . P?j . In fact, it was anticipated by the physicists (see and m the context 
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of Mayers' theory of cluster integrals as we will see below. The form given here, in the 
structural language of species, is the most general one since all the series expansions follow. 
It is also the easiest form to prove. 

Recall that for any species F = F(X), the derivative F' of F is the species defined as 
follows: an F'-structure on a set U is an F-structure on the set U U {*}, where * is an 
external (unlabelled) element. In other words, one sets 

F'[U] =F[U + {*}]. 

Moreover, the operation F i— ► F' , of pointing (or rooting) .F-structures at an element of 
the underlying set, can be defined by 

F* = X ■ F'. 

Theorem 1.1 Let B be a class of 2-connected graphs and C B , the species of connected 
graphs all of whose blocks are in B. We then have the functional equation 

C' B = E{B'(C' B )). (7) 




Figure 4: C' B = E{B'{C B )) 

Proof. See Figure EJ □ 
Multiplying ((JJ) by X, one finds 

C B = X-E(B'(C B )), (8) 

and, for the exponential generating function, 

C' B (x)=x-eMB'(C B (x))). (9) 

1.2 Weighted versions 

Weighted versions of these equations are needed in the applications. See for example 
Uhlenbeck and Ford [^5] • A weighted species is a species F together with weight functions 
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w = wjj : F[U] — > IK defined on F-structures, which commute with the relabellings. Here 
IK is a commutative ring in which the weights are taken, usually a ring of polynomials 
or formal power series over a field of characteristic zero. We write F = F w to emphasize 
the fact that F is a weighted species with weight function w. The associated generating 
functions are then adapted by replacing set cardinalities \A\ by total weights 

\ A \w = ^( a )- 

The basic operations on species are also adapted to the weighted context, using the concept 
of Cartesian product of weighted sets: Let (A, u) and (B, v) be weighted sets. A weight 
function w is defined on the Cartesian product A x B by 

w(a, 6) = u(a) ■ v(b). 

We then have \A x B\ w = \A\ U ■ \B\ V . 

Definition. A weight function w on the species Q of graphs is said to be multiplicative 
on the connected components if for any graph g G G[U], whose connected components are 
ci,c 2 , . . . ,Cfc, we have 

w (g) = w(cx)w(c 2 ) ■ ■ ■ w{c k ). 

Examples 1.2. The following weight functions w on the species of graphs are multiplica- 
tive on the connected components. 

1. wi(g) := y e ^ 9 \ where e{g) is the number of edges of g. 

2- W2{g) = graph complexity of g := number of maximal spanning forests of g. 

3- ws(g) := x^x^x^ 2 ■ ■ •, where n« is the number of vertices of degree i. 

Theorem 1.2 Let w be a weight function on graphs which is multiplicative on the con- 
nected components. Then we have 

g w = E (c w ) . (io) 

For the exponential generating functions, we have 

G w (x) = exp(C w (x)), 

where G w (x) = £„>o l</MUf£ = E n >o(E ge G[n] w (9))%, an d similarly for C w (x). 

Definition. A weight function on connected graphs is said to be block-multiplicative if 
for any connected graph c, whose blocks are b\ , 62 , ■ ■ ■ , bj~ , we have 

w(c) = 10(61)10(62) • • • w(bk). 

Examples 1.3. The weight functions wi(g) = y e ^ and u>2(<?) = graph complexity of 
g of Examples 1.2 are block-multiplicative, but not the function w^(g) = x^x^x 1 ^ 2 ■ ■ ■. 
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Another example of a block-multiplicative weight function is obtained by introducing 
formal variables yi (i > 2) marking the block sizes. In other terms, if the connected graph 
c has rii blocks of size i, for i = 2, 3, . . ., one sets 



(11) 



The following result is then simply the weighted version of Theorem 1.1. 

Theorem 1.3. Let a; be a block-multiplicative weight function on connected graphs 
whose blocks are in a given species B. Then we have 



1.3 Outline 

In the next section, we see how equations (jl()(l and (|12j) are involved in the thermodynam- 
ical study of imperfect (or non ideal) gases, following Mayers' theory of cluster integrals 
|21j . as presented in Uhlenbeck and Ford |29j . In particular, the virial expansion, which is 
a kind of asymptotic refinement of the perfect gases law, is established rigourously, at least 
in its formal power series form. See equation (|34|) below. It is amazing to realize that the 
coefficients of the virial expansion involve directly the total valuation |B a [n]| TO , for n > 2, 
of 2-connected graphs. An important role in this theory is also played by the enumerative 
formula (|42|) for labelled Husimi trees according to their block-size distribution, which 
extends Cayley's formula n n ~ 2 for the number of labelled trees of size n. 

Motivated by this, we first consider, in Section 3, the labelled enumeration of some 
classes of connected graphs of the form Cg, according or not to their block-size distribution. 
Included are the species of Husimi graphs, cacti, and oriented cacti. The methods involve 
the Lagrange inversion formula and Priifer-type bijections. It is also natural to examine the 
question of unlabelled enumeration of these structures. This is a more difficult problem, 
for two reasons. First, equation (|12[) deals with rooted structures and it is necessary to 
introduce a tool for counting the unrooted ones. Traditionally, this is done by extending 
the Dissimilarity charactistic formula for trees of Otter [22] • See for example |Hj • Inspired 
by formulas of Norman (j^j, (18)) and Robinson (j2Zj> Theorem 7), we have given over 
the years a more structural formula which we call the Dissymmetry theorem for graphs, 
whose proof is remarkably simple and which can esily be adapted to various classes of 



tree-like structures; see j2], |H]) CLU US]- [TH] , [TH] . [TH] , Second, as for trees, it should 



not be expected to obtain simple closed expressions but rather recurrence formulas for the 
number of unlabelled Cg-structures. Three examples are given in that section. 

Finally, in Section 4, we present the molecular expansion of some of these species. This 
expansion can be computed first for the rooted species, by iteration, and the dissymmetry 
theorem is then invoked for the unrooted ones. The computations can be carried out 
easily using the Maple package "Devmol" developped at LaCIM and available at the URL 
www.lacim.uqam.ca; see 

Acknowledgements. This talk is partly taken from my student Melanie Nadeau's 
" Memoire de maitrise" |23j . I would like to thank her and Pierre Auger for their consider- 
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able help, and also Abdelmalek Abdesselam, Andre Joyal, Gilbert Labelle, Bob Robinson, 
and Alan Sokal, for useful discussions. 



2 Some statistical mechanics 



2.1 Partition functions for the non-ideal gas 

Consider a non-ideal gas, formed of N particles interacting in a vessel V C IR 3 (whose 
volume is also denoted by V) and whose positions are xt, x~2, . . . , x]}. The Hamiltonian of 
the system is of the form 

+ u(xty\+ Yl v(l^-^l). (13) 

/ l<i<j<N 



H 




where pi is the linear momentum vector and is the kinetic energy of the i th particle, 
U(x~i) is the potential at position xt due to outside forces (e.g., walls), \x\ — Xj\ = is 
the distance between the particles x~t and Xj, and it is assumed that the particles interact 
only pairwise through the central potential <p(r). This potential function ip has a typical 
form shown in Figure a). 





Figure 5: a) the function <p(r), b) the function /(r) 

The canonical partition function Z(V,N, T) is defined by 

Z(Y, N, T) = J exp (-(3H)dF, (14) 

where h is Planck's constant, /3 = , T is the absolute temperature and k is Boltzmann's 
constant, and T represents the state space x±, . . . , xjy, pi, . . . ,p~jy of dimension 6iV. A first 
simplification comes from the assumption that the potential energy U{x\) is negligible or 
null. Secondly, the integral over the momenta pf in (|14|) is a product of Gaussian integrals 
which are easily evaluated so that the canonical partition function can now be written as 

Z(V,N,T) = jyi^Af J ■ " J v exp I -^XJ^d^i - l^i ' --dxN, ( 15 ) 
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where A = h(2irmkT)~2. 

Mathematically, the grand- canonical distribution is simply the generating function for 
the canonical partition functions, defined by 

oo 

Z gI (V,T,z) = Z(V,N,T)(\ 3 z) N , (16) 

where the variable z is called the fugacity or the activity. All the macroscopic parameters 
of the system are then defined in terms of this grand canonical ensemble. For example, 
the pressure, P, the average number of particles, N, and the density, p, are defined by 



I- = ±-\ogZ gT (V,T,z), N = z^- log Z gI (V,T,z), and p := y 



2.2 The virial expansion 

In order to better explain the thermodynamic behaviour of non ideal gases, Kamerlingh 
Onnes proposed, in 1901, a series expansion of the form 

P N , s (N\ 2 , . /7V\ 3 , . 

kf=v + ^ T) {v) + ^ T) {v) + -"' (18) 

called the virial expansion. Here "f2(T) is the second virial coefficient, 73 (T) the third, etc. 
This expansion was first derived theoretically from the partition function Z gT by Mayer 
|21j around 1930. It is the starting point of Mayer's theory of "cluster integrals". Mayer's 
idea consists in setting 

l + f ij =exp(-p^(\xt-x}\)), (19) 

where = f(rij). The general form of the function /(r) = exp(— /3<p(r)) — 1 is shown 
in Figure b). In particular, f(r) vanishes when r is greater than the range n of the 
interaction potential. By substituting in the canonical partition function (|15fl . one obtains 

Z(V,N,T) = —L W [ ■■■ f ]J {l + fa)<m---m. (20) 

JV JV l<i<j<N 

The terms obtained by expanding the product rii<i<j<Ar(l + fij) can be represented 
by simple graphs where the vertices are the particles and the edges are the chosen factors 
fij. The partition function (|20|) can then be rewritten in the form 



dx~N 



9 mw {ij}6<; 



g&5[N) 



E W{g), (21) 



where the weight W(g) of a graph g is given by the integral 

W{g) =/ ■■■/ II AjdSj ■■■dm- (22) 



{i,j}&9 
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For the grand canonical function, we then have 



Z p (V,T,z) = JTZ(V,N,T)(\ 



3 z) N 



N=0 
oo 



< - 3\N 



E mm E w{g){\* z ) 

N=o geQ[N] 
N=0 ' g£Q[N] 

Qw(z). (23) 



z N 



Proposition 2.1 The weight function W is multiplicative on the connected components. 
For example, for the graph g of Figure we have 

W(g) = / fvifnf-nhhfwfmfsndxl ■■■dxs 

f \2 fn f2idxldx2dxj J dx% J U5 f 46 f 56 fssd^dxsdxedxs 
= W( Cl )W(c 2 )W(c 3 ), 

where c\, c 2 and C3 represent the three connected components of g. Following Theorem 
1.2, we deduce that 

Qw{z) = exp(C w (z)), (24) 
where Cw denotes the weighted species of connected graphs, with 

z n 

C w (z) = Y,\C[n]\w— v 



n>l 



and 



\C[n]\ W = E Sy'"j v II fij dxt---dx£. (25) 



ceCH " {i,j}ec 



Historically, the quantities b n (V) = \C[n] \w are precisely the cluster integrals of Mayer. 
Equation (|24|) then provides a combinatorial interpretation for the quantity ^f. Indeed, 
one has, by 1)17(1 . 



P 1 

kf V 
1 



log Z gI (V,T,z) 



log Gw(z) 

^C w (z). (26) 
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Proposition 2.2 For V large, the weight function w{c) = yW(c), defined on the species 
of connected graphs, is block-multiplicative. 

Proof. First observe that for any connected graph c on the set of vertices [k] = {1,2,... k}, 
the value of the partial integral 

/ = I(xt) = Km /"•••/ J] fijdxl ■ ■ ■ dxj~l (27) 
v->oo J v Jv r ..r_ 

is in fact independent of x%. Indeed, since the /^'s only depend on the relative positions 
fij = \xt — Xj\, and considering the short range r\ of the interaction potential and the 
connectednes of c, we see that the support of the integrand in (|27j) lies in a ball of radius 
at most [k — l)r\ centered at and that a simple translation aJf i— >> xf + it will give the 
same value of the integral. It follows that 



W(c) = ■■ fijdx{ ■ ■ ■ dxJZ{ / dxt 

JV JV r f c Jv 



{«j'}Gc 

for V large, and the value of the partial integral (|27|) is in fact w(c): 

lim / ••• / fijdxt ■ ■ ■ dxj^t = lim -W(c)=w(c). (28) 

{«j}ec 

Now if a connected graph is decomposed into blocks &i, &2, • • • , we have 

u;(c) = w{b\)w{b2) ■ ■ ■ w{bk)- 
For example, for the graph c shown in Figure El we have 

w ( c ) = / f 12 f 13 f 23 fzifttf 37 f^f^f 6% f7%dx{dx2 • • • dxj 

Jv 7 

= J fnfvif-adxldx^ fzidxl f 56 dx& f 37 f 36 f 67 f 68 frsd^dxtdx^ 
= w(bi)w(b2)w(b3)w(b4). 

□ 

By Theorem 1.3, it follows that 

C' W = X-E(B' W (C'J), (29) 

where B = B a is the species of all 2-connected graphs, and for the exponential generating 
functions, 

C' w (z) = zeMB' w (C' w (z))). (30) 
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Figure 6: A connected graph with blocks 62, 63, 64 



2.3 Computation of the virial expansion 

The virial expansion (|18|) can now be established, following Uhlenbeck and Ford |29| , From 
(flTf) . we have, for the density p{z) 



N 

v ■ 



d 1 

= z-^-C w (z) 
az 

= CM- (31) 
Hence p(z) satisfies the functional equation (|30|) . that is 

p(z)=zexpBMz))- (32) 



We have, by © and J3TJ) 

P 1 



The idea is then to use this relation in order to express z in terms of p in -^(z), as follows. 



log Z gr (V,T,z 
C w (z) 

C' w (t) dt 
Pit) 



t 



dt. (33) 



Let us make the change of variable 

t = t(r) =rexp(-B' w (r)), 

which is the inverse function of r = p(t), by (|32|) . Note that p(0) = and p(z) = p, and 
also that 

dt = [exp (-£») - r exp • fi£(r)]dr. 
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Pursuing the computation of the integral (|33|) , we have 



p - rm M 



kT Jo t 

rp 



1 — rB'^(r)) dr 
rp 

p— rB'w(r) dr 
Jo 

J ° n>l H - 

n 

p-^(n-l)^, (34) 



ni 

n>2 



N 
V 



where we have set B w (r) = J2 n >2 Pn^r- This is precisely the virial expansion, with p 
Hence the n th virial coefficient, for n > 2, is given by 

IniT) = - (J n 

n! 

= -^^\B[n]\ w . (35) 

Mayer's original proof of the virial expansion is more technical, since he is not aware 
of a direct combinatorial proof of equation ()3U|) . The following observation is used: By 
grouping the connected graphs on [n] whose block decomposition determines the same 
Husimi graph on [n], and then collecting all Husimi graphs having the same block-size 
distribution, one obtains, using the 2-multiplicativity of w, 

\C[n]\ w = w ( c ) 

ceC[n] 

e e n«fo) 

{Bi,B2,.-}eHu[n] {bieB[Bi]} i 

e n e 

{Bi,B 2 ,...}€Hu[n] i b£B[Bi] 

e n^i 

{Bi,B 2 ,...}6Hu[n] i 

]T HuKn 3) ...)^3 3 -- 

n 2 ,n 3 ,... 

^ n i (i-l)=n-l 

where Hu[n] denotes the set of Husimi graphs on [n] and Hu(ri2,n3, . . .) is the number of 
Husimi graphs on [n] having m blocks of size i, for % > 2. Mayer then proves the following 
enumerative formula 

, , (n-ljlnS^" 1 , . 

^'^■••^u^W 1 (36) 
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m 



o 




P 



Figure 7: A Husimi graph with block-size distribution (4, 1, 2, 1, 0, . . .) 

and goes on proving <j5U|l and Q34JI analytically. 

Formula (|36j) is quite remarkable. It is an extension of Cayley's formula n n ~ 2 for the 
number of trees on [n] (take ni = n — 1, 713 = 0, . . .). It has many different proofs, using, 
for example, Lagrange inversion or a Priifer correspondence, and gives the motivation for 
the enumerative problems related to Husimi graphs, cacti, and oriented cacti, studied in 
the next section. 

2.4 Gaussian model 

It is interesting, mathematically, to consider a Gaussian model, where 



which corresponds to a soft repulsive potential, at constant temperature. In this case, all 
cluster integrals can be explicitly computed (see [22] ): The weight w(c) of a connected 
graph c, defined by (|28|). has value 



where e(c) is the number of edges of c and 7(c) is the graph complexity of c, that is, 
the number of spanning subtrees of c. This formula incorporates three very descriptive 
weightings on connected graphs, which are multiplicative on 2-connected components, 
namely, 



wi(c) = y e(c) , io 2 (c) = 7(c), 
already seen in Examples 1.2, and 103(0) = ii n_1 , where n is the number of vertices. 




(37) 




(38) 
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3 Enumerative results 



In this section, we investigate the enumeration of Husimi graphs, cacti and oriented cacti, 
according, or not, to their block size distribution. Recall that these classes can be viewed 
as species of connected graphs of the form Cq and that the functional equation (JSJ) for 
rooted Cg-structures can be invoked, as well as its weighted version (|12|). 

3.1 Labelled enumeration 

Proposition 3.1 The number H n = | Hu[n]| of (labelled) Husimi graphs on [n], forn > 1, 
is given by 

R n = J2S(n-l,k)n k -\ (39) 
fc>0 

where k represents the number of blocks and S(m, k) denotes the Stirling number of the 
second kind. 

Proof. The species Hu of Husimi graphs is of the form Cq with B = K>2, the class 
of complete graphs of size > 2. From the species point of view, it is equivalent to take 
B = -E>2, the species of sets of size > 2. We then have B' = E>\, with exponential 
generating function E>i(x) = e x — 1. Hence the species Hu* of rooted Husimi graphs 
satisfies the functional equation 

Hu* = XS(%(Hu')). (40) 
which, for the generating function Hu*(x), translates into 

Hu*(x) = xR(Ru'(x)), (41) 

with R(x) = exp(e x — 1). The Lagrange inversion formula then gives 

[/]Hu*(x) = i[t n-1 ] (e e * -1 )" 

= I[ t n-l] e n(e*-l) 



n 



n to k] 

fc>o n - 

Since we are dealing with exponential generating functions, we should multiply by n! to 
get the coefficient of . Moreover we should divide by n to obtain the number of unrooted 
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Husimi graphs. In conclusion, we have H n = « n '[^ n ] Hu*(x) = J2k>o S(n — 1, k) n k ~ l . The 
fact that k represents the number of blocks will appear more clearly in the bijective proof 
given below. □ 

We now come to Mayer's enumerative formula 1)361) . 

Proposition 3.2 (Mayer |22j . Husimi ^2j) Let (n2, n^, . . .) be a sequence of non-negative 
integers and n = J2i>2 n «(* — 1) + 1- Then the number Hu(?i2, n%, . . .) of Husimi graphs on 
[n] having raj blocks of size i is given by 

Hu(n 2 ,ra 3 ,...) = — ^~^ ! — : n k -\ (42) 

v ' iJ ' (l!)™2n 2 !(2!)™3n 3 !--- ' v ; 

where k = ~}Zi>2 n i ^ s the total number of blocks. 

Proof. Here we are dealing with the weighted species Hu„, of Husimi graphs weighted by 
the function w(h) = y^V^ • ■ ■ describing the block-size distribution of the Husimi graph 
h. Technically, we should take B w = J2m>2(.E m ) ym , where the index y m indicates that the 
sets of size m have weight y m , for which 



B v{x) = V y m —r and B' v {x) = V y m +i : 

m - ^1 

m>2 m>l 

The functional equation (|12[) then gives 

x m 



Hu^(x) =x(exp(^ y m+1 —) oRu' w (x)) (43) 



Jm+i — : 
ml 

m>l 



and Lagrange inversion formula can be used. We find 



n 

1 r.r,_1i I , ^ > t 



[x n ]Rn' w (x) = _[f-i] exp(5> m+1 — ) 



n \ "„ m\ 

\ m>l 

1 t m 
-[t n-1 ]exp(n V y m +i— r) 
n ^rr, m! 

m>l 

n f-f k\ I ^— ' m! / 

fe>0 \m>l / 



A-2 



1 __ n fc / f i 2 t 3 V 
^ \ETr K + ^2! +y4 3!--j 

K>0 \ / 

>-^t k £ t ^(4)"(4) 

fc>0 fci+fe 2 +— =fe \ / 

n £? n , , , h\k 2 \--- (l!) fe i(2!) fc 2 ••• 

fc>0 k 1 +k 2 +---=k 1 z \ / \ J 
\ ~« \ n ^ fei fc 2 

fc 1 +2fc 2 H =n— 1 
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Extracting the coefficient of the monomial y% 2 y% 3 • • • and multiplying by (n — 1)! = -n! 
then gives the result. □ 

There are many other proofs of (|42|) . Husimi ^2] establishes a recurrence formula 
which, in fact, characterizes the numbers Hu(ra2, n^, . . .). He then goes on to prove the 
functional equation (|43|) . Mayer gives a direct proof which becomes more convincing when 
coupled with a Priifer-type correspondence. Such a correspondence is given by Springer in 
|28j for the number Oc(ra2, n^, . . .) of labelled oriented cacti having block size distribution 

It is easy to adapt Springer's bijection to Husimi graphs. To each Husimi graph h on 
[n] having k blocks is assigned a pair (A, 7r), where A is a sequence . . . ,jk-i) of elements 
of [n] of length k — 1, and tt is a partition of the set [n]\{jfc_i} into k parts. Moreover, if 
h has block size distribution 2 n2 3™ 3 • • •, then tt has part-size distribution l n2 2™ 3 • • •. This 
is done as follows. A leaf-block b of a Husimi graph h is a block of h containing exactly 
one articulation point, denoted by %■ Let b = b{h) be the leaf-block of h for which the set 
b{h)\{3b(h)} contains the smallest element among all sets of the form b\{jb}. 

The correspondence proceeds recursively by 

1. adding to the sequence A, 

2. adding the set b(h)\{jf,^} to the partition tt, 

3. removing the block b(h) (but not the articulation point ju^) from h, and 

4. pursuing with the remaining Husimi graph. 

The procedure stops after the (k — l) th iteration when, in supplement, the last remaining 
block minus the (k — l) th articulation point jk-i is added to the partition n. An example 
is shown in Figure |H1 This procedure can easily be reversed and the resulting bijection 
proves both (JHSJ) and (g^J). 

We now turn to the species of oriented cacti, defined in Example 1.1.6. These struc- 
tures were introduced by Springer in [2B] for the purpose of enumerating "ordered short 
factorizations" of a circular permutation p of length n into circular permutations, Oj of 
length i for each i. Such a factorization is called short if J2i>2(^ ~ l) a « = n ~ 1- 

Proposition 3.3 The number Oc n = | Oc[n]| of oriented cacti on \n\, for n > 2, is given 
by 

where k = X)i>2 n i ^ s the number of cycles. 

Proof. The proof is similar to that of Q39J) . Here B = J2k>2^k, the species of oriented 
cycles of size k > 2 and B' = J2k>iX k , the species of "lists" (totally ordered sets) of 
size k > 1. One can use Lagrange inversion formula, with R(x) = exp(y^). Alternately, 

observe that the factor Ck-i) ™ @5t represents the number of partitions of a set of 
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(3,8) 

{{1.6,9,11}, {2,7}} 





(3,8,8) 
{{1,6,9,11}, {2,7},{3}} 



(3,8,8,4) 
{{1,6,9,11},{2,7},{3},{5,10,12}} 



(3,8,8,4) 

{{ 1,6,9,1 1},{2,7},{3}, {5, 10,12}, {8, 13}} 



Figure 8: Priifer correspondence for a Husimi graph 



size n — 1 into k totally ordered parts so that the Priifer- type bijection of Springer [2S] 
can be used here. □ 



Proposition 3.4 [2S] Let (77,2,^3,...) be a sequence of non-negative integers and n = 
J2i>2 n i{'i — 1) + 1- Then the number Oc(n 2 , n^, . . .) of oriented cacti on [n] having rii 
cycles of size i for each i, is given by 

Oc(n 2 ,n 3 ,--0= (r T, 1)! ^ ( 45 ) 
n 2 !n 3 ! • • • 

where k = J2i>2 n i * s the number of cycles. 

Proof. Again, it is possible to use Lagrange inversion or the Priifer-type correspondence 
of Springer. However the result now follows simply from equation 1)42 J) since it is easy to 
see that 

Oc(n 2 , n 3) • • •) = EK* - ^ Hu ( n2 ' n3 ' • • ■)■ 

i>2 

Indeed, a set of size i can be structured into an oriented cycle in (i — 1)! ways. □ 
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Finally, let us consider the species Ca of cacti which is of the form where B = 
J2k>2 Pk i s the species of polygons. By convention, a polygon of size 2 is simply an 
edge (K2 = E2). See Example 1.1.3 and Figure El These structures frequently appear 
in mathematical research, for example, more recently, in the context of the Traveling 
Salesman Problem (see for instance [3]) and in the classification of graphic matroids [?]. 

Proposition 3.5 (Ford and Uhlenbeck Let (712,71.3, . . .) be a sequence of non-negative 
integers and n = J2i>2 n «(* — 1) + 1- Then the number Ca(n2, 713, . . .) of cacti on [n] having 
Hi polygons of size i for each i, is given by 

Ca(n 2 , n 3 , • • •) = "v^ W^l ^ ( 46 ) 

where k = X)i>2 n i ^ s ^ e number of polygons. 

Proof. Indeed, since any (labelled) polygon of size > 3 has 2 orientations, we see that 

Oc(n 2 , n 3 ...)= 2^>3 n J Ca(n 2 , n 3 , • ■ ■)■ 
The result then follows from (|45|) . □ 

As observed in [£], this corrects the formula given in the introduction of which 
is rather Husimi's formula (|42|) for Husimi graphs. By summing over the polygon-size 
distribution, we finally obtain: 

Proposition 3.6 The number Ca n = | Ca[ra]| of (labelled) cacti on [n], for n > 2, is given 
by ^ ^ 

fc>0 n 2 +n 3 +-=k z n 2 -ri3- 

Tig + 2 7^3 H — n — 1 

3.2 Unlabelled enumeration 

For unlabelled enumeration, the cases of rooted and of unrooted Cg-graphs must be treated 
separately. There is a basic species relationship which permits the expression of the 
unrooted species in terms of the rooted ones. It plays the role of the classical Dissimilarity 
characteristic theorem for graphs (see ^Oj). Let B be a species of 2-connected graphs and 
C73, the associated species of C/j-graphs, that is of connected graphs with blocks in B. We 
introduce the following notations: 

1. C% is the species of Cg-graphs with a distinguished block, 

2. Cf is the species of Cg-graphs with a distinguished vertex-rooted block. 

Theorem 3.7 Dissymmetry Theorem for Graphs ^H]; jl^j - [2| The species C& of con- 
nected graphs whose blocks are in B and its associated rooted species are related by the 
following isomophism: 

C'b + C% = C b + C$. (48) 
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This identity can also be written as 

C B + B{C B ) = C B + C' B -B'(C B ). (49) 

Proof. The proof of (|48[) is remarkably simple. It uses the concept of center of a Cs- 
graph, which is defined as the center of the associated block-cutpoint tree (see Figure 12 
c)). The center will necessarily be either a vertex (in fact an articulation point), or a block 
of the Cg-graph. Now a structure s belonging to the left-hand-side of (|1S|) is a C#-graph 
which is rooted at either a vertex or a block (a cell). It can happen that the rooting 
is performed right at the center. This is canonically equivalent to doing nothing and is 
represented by the first term in the right-hand-side of (|4H|). On the other hand, if the 
rooting is done at an off-center cell, a vertex or a block, then there is a unique incident 
cell of the other kind (a block for a vertex and vice-versa) which is located towards the 
center, thus defining a unique C^-structure. It is easily checked that this correspondence 
is bijective and independent of any labelling, giving the desired species isomorphism. 

For (|4U)) . it suffices to verify the species isomophisms C B = B{C B ) and = C B ■ 
B'{C B ). Details are left as an exercise. □ 

The method of enumeration of unlabelled Cg-graphs then consists in enumerating first 
the rooted unlabelled Cg-graphs, using the functional equation (jHJ or (fl"2|) . and then the 
unrooted ones, using ()49j) . Notice that these steps will lead not to explicit but rather to 
recursive formulas for the desired numbers or total weights. Here we illustrate the method 
in detail in the case of triangular cacti (Example 1.1.4). We will also give some results for 
Husimi graphs and oriented cacti. 




Figure 9: Unlabelled triangular cactus 

Recall that a triangular cactus is a connected graph all of whose blocks are triangles. 
See Figure for an example. These structures (and also the quadrangular cacti) were 
enumerated by Harary, Norman, and Uhlenbeck Here we go further by giving 
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recurrence formulas for their numbers. Let 5 = and A = 5* denote the species of 
triangular cacti and of rooted triangular cacti, respectively. Also set 

S(x) = J2 d nX n and A(x) = J2 D nX n , (50) 

n>l n>l 

where d n and D n denote the numbers of unlabelled triangular cacti and rooted triangular 
cacti, respectively. Here we have B = K 3 = E3, and B' = E 2 . The functional equations 
© and (gnj give 

A = XE{E 2 (A)) (51) 

and 

A + E 3 (A) = 5 + AE 2 (A), (52) 

respectively. Note that 

ZeoE 2 (xi,x 2 ,x 3 , . . .) = Z e (Ze 2 (xi,x 2 , ■ ■ ■), Z E2 (x 2 ,X A , ...),...) 

= Z E (-(xl + x 2 ),-(xl + x 4 ),...) 

( ^ 1 {x\ + x 2k ) 

From (|51|). we deduce (see fTT] ) 



A(x) = (XE(E 2 (A))r(x) 

= xZ EoE2 (A(x),A(x 2 ),A(x 3 ),...) 



= xexp^^(A 2 (x fe ) + A(x 2fc ))y (53) 

A recurrence formula for D n can be obtained by taking x x the derivative of (|53j) . First 
set 

K*) = ^ E ^ (a V) + A ( * 2fc ) 

k>l 



We then have xA'(x) = A(x) + A(x)6(x) and 

xA'(x) - A(x) = A{x)b(x). (54) 



But 



(ix 2/c 

fc>i 

1 



fx 2fc ' 



x • V — 2A(x fc )A'(x fc )A;x fc - 1 + A'(x 2fc )2A:x 2fc - 
ifc>i 2fc V 

2fc\ 2fc 



E (A(x fe )A'(x fe )x fc + A'(x 2fc )x 



fc>i 
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fc>l \h>lj>l j>l 



E IEp^+ E 

m>l \ d|m j=l c£|m 

\ d> 2 d even 



By extracting the coefficient of in (j54j) . we obtain 

Proposition 3.8 The numbers D n , of unlabelled rooted triagular cacti on n vertices sat- 
isfy D\ = 1 and the recurrence formula, for n > 1, 

-. n Id j \ 

= - E D n~m + 1 E E^d-i^ + E 2^1 ' (55) 

772=1 \ d\m j=l d\m I 

\ d>2 d even / 

In order to enumerate unlabelled unrooted triangular cacti, we use (|49|) . that is, 

5 = A + £ 3 (A) - A£ 2 (A). 
The passage to the tilde generating functions yields (see [H]) 

S(x) = A(x) + (E 3 (A(x))r - A(x)(E 2 (A(x))r 

= A(x) + Z E3 (A(x), A(x 2 ), . . .) - A(x)Z E2 (A(x), A(x 2 ), . . .) 

= A(x) + -(A 3 (x) + 3A(x)A(x 2 ) + 2A(x 3 )) - -A(x)(A 2 (x) + A(x 2 )) 
6 2 

= A(x) + i(A(x 3 )-A 3 (x)). (56) 
By extracting the coefficient of x n , we finally obtain 

Proposition 3.9 The numbers d n , of unlabelled (unrooted) triangular cacti on n vertices 
satisfy, for n > 1, 

d n = D n + i j X(3|n)Z>« - XI D i D j D k ■ (57) 

We now turn our attention to the species Hu of Husimi graphs. Let us denote by h n 
and H n , the numbers of unlabelled Husimi graphs and rooted Husimi graphs, respectively, 
and set 

h(x) = ^ h n x n = Hu(x) and H(x) = ^ H n x n = Ru'(x). (58) 

n>l n>l 

Here, as seen in Section 3.1, B = E>2 and B' = E>i = E + , and the basic functional 
equation (|4T1|) translates, for the tilde generating function, into (see [23], p. 51) 



H(x)=xe^ 1 - exp I £ 1 
fc>l \ \m>l 



E -^(^ mfc ) - 1 • ( 59 ) 
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We introduce the auxiliary series 



b(x) = £ b n x" = 4- £ \ exp fx) -^00 I - 1 I , (60) 



n>l fc>l \ \m>l 



n>l 



,?>4 



^(x) = ^ Vn x n = xexp ^ T H(x j ) . (61) 



Then, after some computations, similar to those for triangular cacti, we find: 

Proposition 3.10 The numbers H n , of unlabelled rooted Husimi graphs on n vertices can 
be computed by the following recusive scheme: ip\ = 1, H\ = 1, and, for n > 1, 

6 « = EE E lH ivh, (62) 

<f|ra fis=l Z | cf — h+1 
1 " 

H n+ \ = — 2_. H n _k+ibk, (63) 
1 " 

<Pn+l = - E Vn-m+l E ( 64 ) 
m=l d\m 

For unlabelled unrooted Husimi graphs, we use the Dissymmetry Theorem, in the form 
(|49j). which gives 

Hu*+£> 2 (Hu*) = Hu + Hu'-£ + (Hu'), 
E + (Ru') = Hu + Hu'-S + (Hu'), 

and finally 

Hu = (£+(Hu')) • (1 - Hu*). (65) 
For the tilde generating series, we deduce that 

h{x) = (exp f £ \H{x k ) ] - 1 ] ■ (1 - 2T(z)) 



•l fc 



x 

and we obtain the following result. 

Proposition 3.11 The numbers h n , of unlabelled Husimi graphs on n vertices satisfy 

n-1 



h n = ip n+1 - X ip k+1 H n _ k , (66) 



fc=i 

where the numbers H n and ip n are given by Proposition 3.10. 
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The same method can be applied to weighted Cg-graphs where the weight is defined 
by that is, with the variables (1/2, 2/3, • • •) marking the block sizes. We have done the 
computations for the weighted species Oc^, of oriented cacti, where B w = J2m>2(Cm)y m 
and B' w = J] m>1 (L m )y m+1 , where (C m ) ym denotes the species of oriented cycles of length 
m and weight y m , and L m = X m is the species of lists of size m, endowed here with 
the weight y m +i- Introduce the following tilde generating series, whose coefficients are 
polynomials in the variables y = (2/2 , U3 , ■ ■ ■ ) : 



0(2; y) := Oc w (x) = ^2o n (y)x n , 

n>l 

0(x;y):=6? w (x) = ^O n (y)/. 

n>l 

The two basic species functional equations are 

OC = X ■ E I J2 (L m ) ym+1 (Ocl) I , (67) 

\m>l I 

C„,(Oc') = Oc+L> 2 ^(Oc'), (68) 
where C w = T, m >l( C m)y m , with C\ = X, y 1 = 1, and L> 2 ,w = E m >2p^ m )ym- Note tliat 

zc w = E ^EM*? = E E y -£m& 

m>l " l d\m d>lh>l UU 

where tfi is Euler's totient function, and that Zl >2 w = J2m>2 VmX™ ■ Prom (|67|). we deduce, 
using the plethystic composition rule for the tilde generating function of weighted species 
(see (4.3.1) of 0), 

0(x; y) = xe X p(E^E Vj+i^ W X*) I , (69) 
\k>i K j>i J 

where y k = (y% , yf, . . .). We then set 

0*(s;y) = (0(s;y))' = E ^(y^, (70) 

n>l 

^;y) = E^(y)x m = 4E7Eii° J ^y fc )- ( 71 ) 

m>l ° X k>l j>l 

From (|68|). we deduce 

o(x; y) = E E ^0(x d ; y d ) - E VmO(x; y) m (72) 
d>l /i>l m>2 

and we obtain, after some computations, the following result. 
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Proposition 3.12 The generating polynomials O n {y) and o n (y) for unlabelled rooted and 
unrooted (resp.) oriented cacti on n vertices are given by the following recursive scheme: 



I n 

O n+ i(y) = - V O n - m+1 (y)b m (y) 

II * — * 



(73) 



m=l 



where 



Mv) = EE E v 



(74) 



d|mi=l 3=1 



and 



n 




(75) 



4 Molecular expansions 

The molecular expansion of a species F is a description and a classification of the unlabelled 
F-structures according to their stabilizers, within the language of species. It is interesting 
and useful to have at hand the first few terms of the molecular expansion of Cg-graphs. 
This is possible with the Maple package Devmol, developped at LaCIM. See [Q. 

As an example, consider again the species Hu w of Husimi graphs, weighted according 
to their block-size distribution. Here B w = E>2w = J2k>2i-^k)y k - This is the species BB 
in the following extract from a Maple session, using the package Devmol. The Devmol 
procedure "CBgraphes(BB,n)" produces the molecular expansion of the species of Cbb- 
graphs, up to the given truncation order n. The expansion is then collected by degree 
(vertex number). The weight variables 1/2, VZi ■ ■ ■ iVn are fi rs t declared to Devmol. After- 
wards, the monomial weights in the y^s appear multiplicatively in the expressions. Here 
is the extract: 

> n := 6; 



n := 6 



> ajoutvv(seq(y [k] , k = l..n)); 



2/1, 2/2, 2/3, 2/4, 2/5, 2/e} 



> BB 



sum(y[k]*E[k] (X) , k = 2..n); 



BB := y 2 E 2 (X) + y 3 E 3 (X) + y 4 E 4 (X) + y 5 E 5 (X) + y e E 6 (X) 



> HuW := CBgraphes(BB,n) : 

> af f ichertable (tablephom(HuW) ) ; 
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1 

X 



2 

V2E 2 {X) 
3 

Xy 2 2 E 2 (X)+y 3 E 3 (X) 
4 

y 2 X 2 y 3 E 2 (X) + X y 2 3 E 3 (X) + y 2 3 E 2 (X 2 ) + y 4 E A (X) 

5 

(2/3 y 2 2 + 2/2 4 ) * 3 + X y 2 2 E 2 (X) 2 y 3 + (y 3 y 2 2 + y 2 4 ) X E 2 (X 2 ) 

+Xy 2 4 E 4 (X)+Xy 3 2 E 2 (E 2 (X))+y 5 E 5 (X)+y 2 X 2 y±E 3 (X) 

6 

y 2 5 X 2 E 2 (X 2 ) + (y 2 2 y A + y 3 y A + y 3 y 2 3 ) X E 2 (X) E 3 {X) 
+2/4 E 2 (X) y 2 2 E 2 (X 2 ) + (y 3 2 y 2 + y 2 5 + 3 y 3 y 2 3 ) X A E 2 (X) 
+y 2 5 E 2 (X 3 ) + (y 3 2 y 2 + y 2 5 ) E 2 (X E 2 (X)) + y 2 y 3 2 X 2 E 2 (E 2 (X)) 
+y 3 y 2 3 E 3 {X 2 ) + y 2 y 5 X 2 E 4 (X) + y 6 E e (X) + y 3 y 2 3 X 2 E 2 (X) 2 
+ 2/3 2/2 3 X e + (y 2 5 + y 2 2 y A ) X 3 E 3 (X) + y 2 5 X E 5 (X) 




Figure 10: Husimi graphs of size 5 
As an illustration, the nine terms of size 5 are shown in Figure 1101 Observe that 
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multiciplicities will frequently occur in these expansions. For size 6, for example, the term 
3^3 y\ X A E2{X) corresponds to the 3 similar types of Husimi graphs shown in Figure HT1 




Figure 11: The three unlabelled Husimi graphs of molecular type y\yz ■ X 4 E2(X) 
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